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Topological phases of matter are primarily studied in systems with short-range interactions. In nature, how¬ 
ever, non-relativistic quantum systems often exhibit long-range interactions. Under what conditions topological 
phases survive such interactions, and how they are modified when they do, is largely unknown. By studying 
the symmetry-protected topological phase of an antiferromagnetic spin-1 chain with l/r a interactions, we show 
that two very different outcomes are possible, depending on whether or not the interactions are frustrated. While 
non-frustrated long-range interactions can destroy the topological phase for a < 3, the topological phase sur¬ 
vives frustrated interactions for all a > 0. Our conclusions are based on strikingly consistent results from 
large-scale matrix-product-state simulations and effective-field-theory calculations, and we expect them to hold 
for more general interacting spin systems. The models we study can be naturally realized in trapped-ion quantum 
simulators, opening the prospect for experimental investigation of the issues confronted here. 
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Since the discovery of topological insulators [1-3], there 
has been tremendous interest in exploring various topological 
phases of matter, both theoretically [4, 5] and experimentally 
[6-8]. Topological phases are generally associated with—and 
derive much of their presumed utility from—stability against 
local perturbations. But precisely what constitutes “local” in 
this context is a subtle issue; power-law decaying (1 /r a ) in¬ 
teractions, which are present in many experimental systems, 
do not necessarily qualify [9-11]. Recent theoretical advances 
have begun to elucidate the conditions under which long-range 
interacting systems maintain some degree of locality [12, 13], 
potentially providing some insight into effects of long-range 
interactions on topological phases of matter. And recently, ex¬ 
plicit theoretical evidence of topological order has been found 
in a variety of long-range interacting systems, including dipo¬ 
lar spins [14] or bosons [15], fermions with long-range pairing 
[16] and hopping [17, 18], and electrons with Coulomb inter¬ 
actions [19]. These results notwithstanding, a complete under¬ 
standing of how topological phases respond to the addition of 
long-range interactions is still lacking. 

The stability of topological phases to small local perturba¬ 
tions is intimately connected to the existence of a bulk exci¬ 
tation gap [20, 21], and the introduction of long-range inter¬ 
actions to a short-range Hamiltonian supporting a topological 
phase poses several potential challenges to this connection. 
First, even if the gap remains finite, long-range interactions 
can change the ground state correlation decay from exponen¬ 
tial to power-law [16, 18, 22, 23]. Thus topological phases 
with local interactions are, at the very least, subject to qual¬ 
itative changes in their long-distance correlations. Second, 
the gap can in principle close in the presence of long-range 
interactions, even when they decay fast enough that the total 
interaction energy remains extensive [20, 24]. Third, long- 
range interactions have the ability to change the effective di¬ 
mensionality of the system [25, 26], and thus might change the 


topological properties even if the gap does not close [16, 18]. 
We emphasize that the understanding of these issues is not of 
strictly theoretical interest. Many of the promising experimen¬ 
tal systems for exploring or exploiting topological phases of 
matter, e.g. dipolar molecules [27-29], magnetic [30] or Ry¬ 
dberg atoms [31], trapped ions [32-37], and atoms coupled 
to multi-mode cavities [38], are accurately described as quan¬ 
tum lattice models with power-law decaying interactions. The 
unique controllability and measurement precision afforded by 
these systems hold great promise to improve our understand¬ 
ing of topological phases [39^-2], but first we must reliably 
determine when—despite their long-range interactions—they 
can be expected to harbor the topological phases that have been 
theoretically explored for short-range interacting systems. 

To address these general questions, in this manuscript we 
study a spin-1 chain with antiferromagnetic Heisenberg inter¬ 
actions, which is a paradigmatic model exhibiting a symmetry- 
protected topological (SPT) phase [43, 44]. Specifically, we 
consider two extensions of the short-range version of this 
model by including long-range interactions that decay either 
as J&{r) = l/r a or as J^(r) = (—1 ) r_1 /r a , which could 
be simulated in trapped-ion based experiments for 0 < a < 3 
[45, 46]. Based on a combination of large-scale variational 
matrix-product-state (MPS) simulations and field-theory cal¬ 
culations, we establish and explain a number of important 
and potentially general consequences of long-range interac¬ 
tions. The J'Jdf) interactions are unfrustrated, being anti¬ 
ferromagnetic (ferromagnetic) between spins on the opposite 
(same) sub-lattice. In this case, numerics and field theoretic 
arguments suggest the destruction of the topological phase for 
a < 3, accompanied by a closing of the bulk excitation gap 
and spontaneous breaking of a continuous symmetry in ID, 
consistent with other recent findings on the relevance of long- 
range interactions for a < D + 2 in U-dimensional quantum 
systems [47, 48]. The J a (r) interactions are frustrated, and, 
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Figure 1: (a) Low-lying energy levels of the Haldane chain for even 
L. The entanglement structure of ground states is shown at the 
bottom. The ground states in the total S z = 0,1, 2 subspace are 
named |0), |1), |2) and have energies E 0 ,Ei 1 E 2 . (b-c) The m th 
largest value A m (m = 1, 2 • • • 8) of the ground-state entanglement 
spectrum for H a (b) and H' a (c) using finite-size MPS calculations 
with L = 200. We choose the |1) state to avoid extra entangle¬ 
ment between edge spins. For H' a , the entanglement spectrum for 
1.5 < a < 4 will exhibit a smooth crossover between the a = 1.5 
and a = 4 cases due to the finite system size, but we expect a sharp 
transition at some a c < 3 in the thermodynamic limit. The exact pair 
degeneracies in {A m } are a result of the spatial-inversion symmetry 
protecting the topological phase [44, 49]. 


remarkably, do not close the bulk excitation gap for any a > 0 . 
In addition, two key properties of the SPT phase, a doubly 
degenerate entanglement spectrum [49] and a non-vanishing 
string-ordered correlation [50], are both preserved. However, 
because of the long-range interactions, spin-spin correlations 
and the edge-excitation amplitudes only decay exponentially 
within some intermediate distance scale, after which they de¬ 
cay algebraically. We expect these qualitative changes to be 
quite general, occurring in other long-range interacting sys¬ 
tems in which the topological phase survives. 

Model —We consider a spin-1 chain with either frustrated 
or unfrustrated long-range Heisenberg interactions: 


these excited states belong to a degenerate ground-state sub¬ 
space in the thermodynamic (L oc) limit. The second set 
of excited states all have S = 2, contain spin excitations in 
the bulk of the chain, and have an energy gap ( bulk gap) that 
converges to a finite value when L oc. The entanglement 
structure of the four ground states is close to that of the Affieck- 
Kennedy-Lieb-Tasaki (AKLT) states [61] shown at the bottom 
of Fig. 1(a), where each spin-1 is decomposed into two spin- 
l/ 2 s, pairs of spin-l/ 2 s on neighboring sites form singlets, and 
the system is finally projected back onto the spin-Is. The four 
quasi-degenerate ground states correspond to the four states 
formed by the two unpaired spin- 1 / 2 s at the edge. 

We use variational MPS calculations [62-65], to determine 
the ground-state entanglement structure of H a and H f a in 
Fig. l(b-c). For a > 0 (a > 3), the ground-state entangle¬ 
ment spectrum of H a (H' a ), defined as the eigenvalues of the 
left/right half-chain’s reduced density matrix, is dominated by 
the two largest degenerate eigenvalues Ai = A 2 ~ 0.5. This 
can be understood heuristically as the result of cutting a spin- 
1/2 singlet in the AKLT state, and suggests the survival of the 
topological Haldane phase. For H' a with a < 3, the entangle¬ 
ment spectrum has an entirely different structure, and we will 
study the related ground state properties below. 

Effective field theory. —The low energy physics of the Hal¬ 
dane chain can be understood via field theoretic analysis due to 
Haldane [52] and Affleck [ 66 ]; here we build on their work to 
provide a field theoretic treatment of the long-range interact¬ 
ing model. We begin by decomposing the spin operators into 
staggered and uniform fields, n(2i + |) = (S 2 i — S 2 i+i)/2 
and l(2i + 5 ) = (S 2i + S 2i+ 1 )/2. The intuition behind this 
decomposition is that the classical ground state of both H a 
and H' a is Neel-ordered for any a > 0, with n 2 (x) = 1 
and l(x) = 0. We therefore expect that in the quantum 
ground state n 2 (x) « 1 , while l(x) ~ 0 represents small 
quantum fluctuations in the direction of n(x). Importantly, 
we expect that only long-wavelength fluctuations of n(x) and 
l(x) will be important at low energy. In momentum space, 
we can write H a « f dq{uj(q)\n(q)\ 2 + Q(q)\l(q)\ 2 ) and 
H’ a « f dq(Q(q)\n(q)\ 2 + u(q)\l(q)\ 2 ) [67], with 

oo oo 

u(q) = 2 E cosqr ’ = 2 E cosqr • ^ 

r= 1 r= 1 


H a =£j a (r)Sj-S j+r , H' a =Y J J»S r S j+r . ( 1 ) 

j,r> 0 j,r>0 

With only nearest-neighbor interactions (a oc), H ^ = 
H ^ is usually called the Haldane chain , which has been exten¬ 
sively studied theoretically [51-53], numerically [54-58], and 
experimentally [59, 60]. The low-lying states of the Haldane 
chain are shown in Fig. 1(a) for an open boundary chain with 
even size L. The unique ground state has total spin 5 = 0. The 
first set of excited states has 5 = 1 (h = 1 ), contains spin ex¬ 
citations only near the edge of the chain, and is separated from 
the ground state by an energy gap (edge gap ) that is exponen¬ 
tially small in L and topologically protected. Consequently, 


For any a > 0, c u(q) is analytic at small q and can be ex¬ 
panded as cco E uj 2 q 2 + 0(q A ), whereas fl(q) is non-analytic 
at small q with an expansion f2 0 EEt 2 q 2 + A|g | a_1 + 0(q 4 ). 
The coefficients cco, 2 , ^ 0 , 2 , and A depend on a , but their exact 
values are not important for the following analysis. Physically, 
the analyticity (non-analyticity) of the spectrum arises because 
the long-range interactions interfere destructively (construc¬ 
tively) for the staggered field. Keeping only the lowest non¬ 
trivial order in q for the dispersion of both n(q) and l(q) 
turns out to be sufficient for obtaining qualitatively correct 
behavior of the excitation gap. Therefore, we keep only the 
0 th order term in the dispersion of l(q), and the next-leading 
term in the dispersion of n(q) (for n(q), the 0 th order term 
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only adds a constant to the Hamiltonian due to the constraint 
n 2 (x) = 1). Thus for a > 0 (a > 3) the Hamiltonian H a 
(H' a ) is approximately given by (ignoring the order-unity co¬ 
efficients) H a ~ H' a ~ f dq(q 2 \n(q )\ 2 + \l(q)\ 2 ). When the 
zero-temperature partition function is expressed as a coherent- 
spin-state path integral, the action is quadratic in the field l 
and it can be integrated out [68, 69]. The remaining path in¬ 
tegral over the staggered field n is a 1+1 D 0(3) nonlinear 
sigma model, with Lagrangian density (nonlinear constraint 
n 2 (x) = 1 implied) 

C(x) « -[\dn/dt \ 2 — v 2 s \dn / dx\ 2 ) . (3) 

Here g is an effective (a- and short-distance-cutoff-dependent) 
coupling strength, and the spin-wave velocity v s is also a- 
dependent. This model is gapped and disordered for all g [51]. 

To investigate ground state properties of Eq. (3), we can re¬ 
move the constraint n 2 (x) = 1, while phenomenologically 
introducing a mass gap A a and a renormalized spin-wave ve¬ 
locity v a (the parameters A' a and v' a will be used to describe 
the Lagrangian for H' a ) [57, 58]. Transforming to momentum 
space, we thereby arrive at a free-field Lagrangian density 

C(q)<x\dn/dt\ 2 -(A 2 a + vlq 2 )\n(q)\ 2 . (4) 


This Lagrangian leads to ground-state correlations Cij = 
(SfSj )o [where (• • -) m denotes the expectation value in the 
state | m) defined in Lig. 1(a)] that decays as 

Cij K G 1 )' f 2 K (—l) r ^o(r/U- (5) 

J y/Al + vlq 2 

Here = v a /A a (or £' a = v' a /A' a for H' a ) defines the corre¬ 
lation length, and Kq(x) is a modified Bessel function, which 
behaves as Kq(x) ~ exp(— x) / A Jx for large x. 

Lor a < 3, the non-analytic \q \ a ~ 1 term in H' a dominates 
the dispersion of n(q) at small q, and Eqs. (3-4) are not valid. 
To analyze this case, we write down the renormalization group 
(RG) flow equation for the coupling strength g under the scal¬ 
ing transformation x -+ xe~ l to one-loop order [68, 70], 


dg a — 3 g 2 
dl 2 ^ 47r 


( 6 ) 


Lor a < 3, an unstable fixed point appears at g* = 2tt( 3 — a ), 
and for a bare coupling g < g* the RG flow is towards a weak- 
coupling ordered state at g = 0 [68]. The bare coupling, and 
therefore the value of a at which this phase transition occurs, 
is difficult to determine a priori. But we nevertheless expect 
(and confirm numerically) that for a < a c , with 2 < a c < 
3, the gap will close as the system spontaneously breaks the 
continuous SU(2) symmetry of H' a [48, 71]. 

Comparison with numerics .—Using finite-size MPS calcu¬ 
lations, we have obtained the bulk excitation gap — E\ and 
the correlation length [fitted using Eq. (5)] for both H a and 
H' a . As shown in Lig. 2(a-b), we see consistent results with 
the field theory predictions. Lor H a , the gap remains open 
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Figure 2: (a) Bulk gap A a and ground state correlation length in 
the L —»• 00 limit, obtained by finite-size scaling for 200 < L < 
500. (b) Bulk gap A' a and £' a with L — 100 and L — 300. (c) 
Ground-state string-ordered correlation function Sij for H a and H' a 
with a = 1.5 and L = 300. For various a and 200 < L < 500, 
we consistently find that Sij quickly saturates to a finite value for H a 
at all a > 0, but vanishes at large distance for H' a at a ;$ 3. (d) 
Ground-state spin-spin correlation Cij for a = 0.5 and L — 500. 
This choice of a s= 0.5 is arbitrary, but assists in a clear presentation 
of the coexisting exponential and l/r a+4 power-law decays. 


for all a > 0, and the correlation length decreases together 
with a due to both an increase of the bulk gap, and a decrease 
of the spin-wave velocity (as a result of a weakened Neel or¬ 
der for longer-range interactions). To the contrary, for H' a , 
the gap decreases quickly as the interactions become longer 
ranged, and the correlation length diverges when a decreases 
to around 3, suggesting the disappearance of the topological 
phase at a < 3 [72]. Calculation of the string-ordered corre- 
lation Sij = {SfSj IIi<fe<j(- 1 ) SJ >o_ of both H a and H' a at 
a = 1.5 [Fig. 2(c)] provides further evidence that the topolog¬ 
ical phase survives for H a , but not for H' a , for 0 < a < 3. 

We now analyze the effects of terms beyond leading order in 
q that have been ignored in our field theory treatment. Includ¬ 
ing the higher-order analytic terms, such as the 0(q 4 ) term, 
will result in negligible corrections to the correlation func¬ 
tions that decay in distance faster than Eq. (5) [57]. However, 
even for a > 3, inclusion of the non-analytic 0(|g| a-1 ) term 
will add a power-law tail to the correlation functions, which 
will dominate over Eq. (5) at long distance. In the supplemen¬ 
tal material, we show by a more involved field-theory calcula¬ 
tion that, for H a , Cij decays as i/r a + 4 at large r. Our MPS 
calculations using L = 500 spins [Fig. 2(d)] show remark¬ 
able agreement with the field theory predictions, even captur¬ 
ing the oscillations in \Cij \ occurring at intermediate distance 
where the short-range and long-range contributions to the cor¬ 
relation functions are of comparable magnitude and interfere. 
A power-law tail in Cij should also exist for H' a , but the in¬ 
creased correlation length prevents us from observing its exis- 
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tence clearly for a > 3. 

Edge excited states .—We expect the influence of long-range 
interactions on the edge- and bulk-excited states to be strong 
at small a\ because the topological phase of H' a does not sur¬ 
vive for a < 3, we will focus on H a from now on. Edges can 
be introduced into the field theory by replacing the two end 
spin-l’s with spin-l/2’s, represented by tl (tr) for the left 
(right) edge, resulting in an edge-bulk coupling Hamiltonian 
He = E& 2 1 Si ■ (• r L /(i - 1)“ + t r /(L - i) a ) [57], For the 
edge excited state |1) [Fig. 1(a)], tl,r are polarized in the +z 
direction, and we expect (S z ) to decay away from the ends. 
Solving the free theory defined by Eq. (4) and treating H c us¬ 
ing standard first-order perturbation theory [57], we find that 
(■ n z (x))i oc f dq(exp[iq(L - x)] - exp [iq(x - 1)])/(A^ + 
V W) exp[-(L - x)/£ a ] - exp[—(x - l)/£ a ] for even 
L. In addition, (l z (x))\ contributes a power-law correction 
l/(x — l) a + 2 + 1 j{L — x) a+2 for x far away from both 
ends [73]. Our numerical calculation of (S z (x)) i, shown in 
Fig. 3(a), agrees well with a sum of these two contributions, 
clearly exhibiting an exponential followed by 1 /r a + 2 decay. 

The edge gap \Ei — Eq\ can be obtained by using a path 
integral to integrate out the n field [57], resulting in an ef¬ 
fective edge-edge Hamiltonian oc (—1 ) L exp (— L /£ a ) ■tr • tr . 
This scaling is confirmed, at relatively small L , by the nu¬ 
merical results in Fig. 3(b). However, the numerics also re¬ 
veal that at large L the edge gap receives a long-range correc¬ 
tion given by 1 /L a . This remarkably simple result, including 
the unity prefactor, can be understood as follows. The edge- 
excited states behave differently from the bulk-excited states 
due to correlations between the orientations of T\ and T 2 , 
and therefore (Si • Sj) i — (S t ■ Sj)o is very small unless i 
and j are very close to 0 and L, respectively. Thus we have 
Et-Eo* L-“ £*<,•«& • Sj )i - (Si ■ Sj ) 0 ) = 1/L a , 
where the last equality is a sum-rule following from the total 
spin of the ground (S = 0) and edge excited (5 = 1) states. 

Bulk excited states .—As in the short-range Haldane chain, 
the elementary bulk excitations of H a are spin-1 magnons 
[55-57]. Physically, the magnon represents fluctuations in the 
field n, and, from Eq. (4), these fluctuations have a disper- 
sion relation e a (q) = i/Aj + (x a g) 2 « A a + q 2 v 2 a /( 2A a ) 
(valid at small q). The lowest-energy magnon wave-function 
\F 0 (:r) can be extracted from the numerics using the relation 
14>o (i) | 2 ~ | (Sf) 2 — (S z ) 1 1. The presence of long-range inter¬ 
actions gives the magnon an additional potential energy due to 
the edge-bulk coupling Hamiltonian H c , and ^(x) can be ap¬ 
proximately described by the following Schrodinger equation 
(with Dirichlet boundary condition at x = 1, L) 


v 2 d 2 ^t(x) 
2A a dx 2 


+ - 


(x — l) a (L — x) 


^(x)=£'f/(x). (7) 


The kinetic (potential) energy always scales as 1/L 2 (1 /L a )\ 
therefore, for a > 2 and large L , the potential energy can 
be ignored. The ground-state energy £q « u 2 7r 2 / (2A a L 2 ) 
and probability density |^ 0 (^)| 2 ~ (2/L) sin 2 (7rx/L) are 
then identical to those of a particle in a box, as confirmed 
numerically in Fig. 3(c-d). The relation E<i — E\ « A a + 



Figure 3: (a) Distribution of an edge excitation in state |1) for L = 
500 and a — 2. (b) Edge gap | E\ — Eq \ as a function of the chain size 
L for a = 3. (c) Lowest energy magnon probability density distri¬ 
bution for L = 200 and a = 3.0, 0.5. (d) The finite-size correction 
to the lowest magnon excitation energy [see Eq. (7)]. For a = 3, we 
obtain v a = 2.18 and v a /A a « 4.51, in good agreement with the 
~ 4.55 obtained in Fig. 2. 


v^tt 2 /(2A a L 2 ) allows us to obtain both v a and through 
finite-size scaling [Fig. 2(b)], and we confirm that the corre¬ 
lation length determined by = v a /A a agrees with that 
obtained by fitting Cij using Eq. (5). For a < 2, the poten¬ 
tial energy dominates the kinetic energy for large L , and the 
potential can be approximated as harmonic around x = L/2. 
Thus |T r 0 (x)| 2 resembles a Gaussian [Fig. 3(c)], and a simple 
scaling analysis predicts a width 7 oc L 1-a / 2 . In the large-L 
limit, 14/o0r)| 2 become sharply peaked at x — L/2 and, from 
Eq. (7), we expect the bulk gap to scale as A a + (2/L) a , which 
is clearly observed in Fig. 3(d). Since E 2 — E\ = 2 when 
a = 0, it follows that A a ^o = 1, consistent with Fig. 2(a). 

Outlook. —The stability of the topological Haldane phase 
to 1 /r a interactions for all a > 0 is favorable for trapped- 
ion based experiments, as stronger couplings can be achieved 
for smaller a [36, 37]. Moreover, because the correlation 
length shrinks for longer-range interactions, a relatively small 
number of ions will suffice to suppress finite-size effects. 
Probing the topological phase by measuring both Cij and 
Sij with single-site resolution is nearly impossible in typical 
condensed-matter systems, but is quite straightforward in ion 
traps [74]. Based on the generality of our field theory anal¬ 
ysis, we speculate that for generic lattice models, the tails in 
the power-law interactions can possibly destroy the topologi¬ 
cal phase only when long-range interactions are unfrustrated 
and a < D + 2. Experimentally, unfrustrated long-range in¬ 
teractions can be easily implemented by generating a 1 /r a fer¬ 
romagnetic interaction [71]. We hope that our work can serve 
as a springboard for future studies on how distinct topological 
phases behave in the presence of long-range interactions. 
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